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TECHNICAL NOTE 2512 


STRESSES IN A TWO -BAY NONCIRCULAR CYLINDER 
UNDER TRANSVERSE LOADS 
By George E. Griffith 


SUMMARY 


A method, taking into account the effects of flexihllity and based 
on a general eighth-order differential equation, is presented for finding 
the stresses in a two-hay, noncircular cylinder the cross se^lon of ' 
which can he composed of circxilar arcs. Numerical examples are given 
for two cases of ring flexihllity for a cylinder of douhly symmetrical 
(essentially elliptic) cross section, subjected to concentrated radial, 
moment, and tangential loads. The results parallel those already obtained 
for shells with circular rings. 


INTRODUCTION 


In airplane fiiselages with flexible rings subjected to concentrated 
loads, the stresses in the neighborhood of the load differ markedly from 
those given by the simple engineering formulas, and more refined methods, 
which take into account the interaction of rings and shell, are rreeded 
to predict the stresses accurately. The first paper on this subject, 
published in May 19^^^ "was that of Wignot, Combs, and Ensrud (reference 1), 
who treated the circular cylinder subjected to concentrated loads but 
neglected the effect of the extensional deformations of the shell. 

Hoff (reference 2) gave a more ccanplete analysis. Including the effects 
of many rings, for the case of symmetric transverse loads. The results 
were corroborated experimentally by Kuhn, Duberg, and Griffith (refer- 
ence 3)^ who also extended the theory to include concentrated moment 
and tangential loads. Later, Duberg and Kempner (references L and 5) 
reduced the labor of computation by giving the results in the form of 
charts and showed that for practical purposes it was usually sufficient 
to consider only a region within 2 bay lengths of the load. F\nrther 
Investigations considered additional effects, heretofore neglected, such 
as the shearing and axial deformation of the rings (references 6 to 8), 
bending rigidity of stringers (references 8 and 9)} shear carried by the 
stringers (reference J), and eccentricity of ring and sheet (references 6, 
7, 8, and 10). 
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All investigations referred .to dealt exclusively with reinforced 
monocoque cylinders of circular cross section. The present paper gives 
an analysis for a two-hay noncircular cylinder, enclosed between ring 
bulkheads rigid in their planes, with the middle, flexible ring subjected 
to concentrated and distributed loads. The fundamental assumptions used 
in the analysis are the same as those previously used for clrcxilar 
cylinders. 

Many noncircular fuselages can be closely approximated by using 
circular sections of different radii and Joining the sections at~ points 
of tangency. The rings discussed herein are of this form. Associated 
with each ring section is a two-bay panel (fig. 1 ), any number- of similar 
panels ccmprising the complete two-bay cylinder. The solution for the 
stresses in such a structure is based on the development of a general 
eighth-order differential equation, written in terms of the moment at 
the skin center line. A separate differential equation of the same 
general form applies to each curved panel in the structure. Application 
of the correct boundary conditions results in sets of simultaneous 
equations which yield the unknown constants in the moment expression. 

Although a two-bay cylinder does not conform to the usual fuselage 
structure, it is believed, on the basis of compgirisons with some of the 
work previously mentioned, that-the results obtained are indicative of 
those found in more complicated structiu*es. In accordance with the 
findings of prior investigations, shear and axial deformations of the 
loaded ring are neglected, but eccentricity of -ring and sheet is included. 

The nvunerical examples deal with a doubly symmetrical, two-bay 
cylinder of nearly elliptic cross section (fig. 2 ), subjected. to concen- 
trated loads. 


SYMBOLS 


A = 


B = 


t'R^ 

Il3 

Et'R^ 

GtL^ 


Cq coefficient of shear flow in sheet 

Cm 

Cy 


coefficient of ring axial force 
coefficient of ring bending moment 
coefficient of ring shear force 
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E 

G 

H 

I 

L 

M 

Me 


Young's modiilus in sheet and ring, pounds per square inch 

shear modulus in sheet, pounds per square inch 

axial force acting on sheet, pounds 

moment of inertia of ring section, inches^ 

length of hay, inches 

bending moment at sheet center line, inch-pounds 
applied concentrated moment, inch-pounds 


Mr 

P 

R 

Rl, B-2 

T 

V 

a ± hi 


Cl + dii, 

C2 ± d2i 
e 


f 

h 


bending moment about neutral axis of ring, inch-pounds 

applied concentrated radial load, pounds 

radius to sheet center line of circular panel. Inches 

radii of circular panels comprising cylinder, inches 

applied concentrated tangential load, pounds 

shear force, pounds 

complex roots of aiixiliary algebraic equation (used when real 
roots also occur) 

complex roots of auxiliary algebraic equation (used when no 
real roots occur) 

eccentricity of ring and sheet (distance between sheet center 
line and neutral axis of ring) , inches 

distributed radial load acting on sheet, pounds per inch 

distributed axial load acting on sheet, pounds per inch 


1 = 

kj_, kg real roots of auxiliary algebraic equation 


m distributed moment acting on sheet, inch-pounds per inch 

p normal force in stringers, pounds per inch 



k 


NACA TN 2512 


q shear flow in sheet, pounds per inch 

s distance in circumferential direction, inches 

t thickness of all material carrying bending stresses in panel 

if uniformly spread around circumference, inches 

t thickness of sheet, inches 

u displacement of sheet in axial direction, inches 

V displacement of sheet in circumferential direction. Inches 

Vq displacement of sheet in circumferential direction at ring, 

inches 

w displacement of sheet or ring in radial direction, inches 

X distance in axial direction measiu*ed from loaded ring, inches 

7 shear strain in sheet 

Sg normal strain in sheet in circumferential direction 

normal strain in sheet in x-direction 

0 angular distance, degrees or radians 

a stringer normal stress, pounds per square inch 

Cg normal stress in sheet in circumferential direction, pounds 

per 8 quare inch 

T shear stress in sheet, pounds per square inch 


GENERAL DIFFERENTIAL EQUATION 


The basic element of the present analysis is a two-bay panel, as 
shown in figure 1, composed of sheet, longitudinal stiffeners or stringers, 
and transverse stiffening ring sections of constant radius. Distributed 
or concentrated loads are applied in the plane of the middle or flexible 
ring. By Joining several panels at points of tangency, many two-bay 
cylinders of various cross-sectional shapes can be achieved. 

In actual practice the sheet covering is outside the rings, and 
this eccentricity of sheet and ring is henceforth taken into account. 
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For convenience the moment considered is the moment which exists at the 
sheet center line rather than at the neutral axis of the ring. From 
this moment the bending moment in the ring is easily found. 

The underlying assumptions used in the analysis are as follows; 

(1) The shear stress, carried by the sheet alone, may vary in the 
circumferential direction but remains constant in the axial direction. 

(2) The material in the cross section of the panel (sheet and 
stringers) capable of canning normal stresses due to bending of the 
panel is assumed spread around the circumference in _a fictitious sheet 
of thickness t* . 

( 3 ) The loaded ring has no torsional stiffness or bending stiffness 
out of its plane. 

( 4 ) The end ring supports are restrained from deforming in their 

planes but are free to warp out of their planes. 


Under these ass\miptlons, for any panel with constant geometrical 
properties, a general differential equation is developed for the moment 
at the sheet center line. All forces, stresses, and displacements in 
the panel may be obtained from that moment. The general differential 
equation (see the appendix for development) is 


d®M ^ d^ ^ L 6e^A 

deo de^ V / de^ 


A d^ _ 


B - 


2eA 


2 - I- + ^1 


^ + 6a/i - e + = 




F(e) 


(1) 


where 
F(0) = R^ 


dh 


^ + ( 1 - 3 B)^ - 3 B 
cie5 cie-J 



+ (2 - 3B)^ . 6 b ^ ^ - 


Sca/- e e^j dmL 
E\^~R'^^d 0 [ 
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When the ratio e/R, a measiire of the eccentricity of the ring and 
sheet, is very small it can he neglected; equation (1) then 

contains only the two nondimensional parameters A and B, which are 
related to the geometry and physical properties of the panel: A = ■ ■■■■ 

,2 Il3 

and B = The ratio a/b is an index of the over-all flexibility 

GtL^ 


of the structure. When this ratio is very small the ring section is 
rather rigid. A large value of a/b (500^ for example) indicates that 
the ring section is somewhat flexible, causing radical departures of the 
stress distributions from the elementary veilues. In actual practice 
the parameter B varies over a small range (between about 10 and 80 ), 

whereas A is usually much larger and may be as great as 2 x 10*^, or 
even more; consequently, flexibility may be thought of in terms of 
parameter A alone. Values of A less than about 200 usually Indicate 
relatively inflexible rings, and Increasing values indicate increasing 
ring flexibility. 


The right-hand side of equation (l), F(0), contains terms resulting 

from the application of distributed loads; if only concentrated loads 
are present P(0) becomes zero. 


SOLUTION OF GENERAL DIFFERENTIAL EQUATION 


The solution of equation (l) depends upon the nature of the roots 
of the auxiliary algebraic equation associated with the differential 
equation. Of the required eight roots of the algebraic equation, four 
may be real and four caaplex, or all eight may be complex. . (For the 
special case of no eccentricity of ring and sheet, sets of roots for 
several combinations of parameters A and B are given in table I . ) 

If there are foxrr real roots ik 2 and four complex roots 

±(a ± hi), the solution for the bending moment is 


M = 


C^e 


kn0 


-k,0 

+ Cj-je 


Cve 


(a+bl)0 


+ Cyje 


^VIII® 


■(a-bl)e 


+ Cjiie^ + Cpye"^® + 

(a-bl)e ^ p ^-(a+bi)0 
+ ^II® 


+ 


+ Particular solution 


( 2 ) 
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which can be written In real form as 

M = C]_ slnh + C2 cosh alnh cosh k20 + 

slnh aS cos bS + Cg slnh a0 sin b© + Cy cosh a0 sin b© + 

C3 cosh a© cos b© + gj(©) ( 3 ) 

If there are fotir pairs of complex roots +(c]_ + and 

±(c2 ± < 32 ^)# bending moment is 

(ci+dil)© (ci-dil)e -(c^+dii)© -(ci-d3_i)0 

M = C^e + Cj^e + Cjj^e + C^y-e + 

Cyjjje 4. Particular solution (i-) 

which, written in real form, becomes 

M = sinh Cj^© cos d]_© +■ C2 slnh C]_© sin d]_© + cosh C]_6 sin d2_© + 

C4 cosh C3_© cos d^^S + slnh C2® cos 12© + Cg sinh c^B sin d2© + 

Cj cosh C20 sin 6^9 + Cg cosh CgS cos d2© + JJ(®) ( 5 ) 

Since the particular solutions gJ(©) and JJ(©) depend upon the form 

of F(©), no formal solutions are given here. Further remarks concerning 
the solution of equation (l) are confined to the solution in real form, 
either equation (3) or equation (5)* 

Equation ( 3 ) or equation ( 5 ) expresses the moment in any particular 
panel where A, B, and R are constant j such a solution exists for each 
panel constituting the structure to be analyzed. For the case of a 
circular cylinder, solution (3) or (5) is the same as the energy solutions 
of references 2 and 3 for a similar two-bay structvire. 


RELATION BETWEEN THE MOMENT M AND OTHER FORCES AND DISPLACEMENTS 


All forces and displacements in the panel may be exp>ressed in terms 
of the moment given by equation (3) or ( 5 ) and derivatives of this moment. 
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As previously noted, this moment Is the moment at the sheet center line, 
and with the exception of the ring bending moment all other forces and 
displacements are those at the center line of the sheet. When no 
eccentricity is Involved the sheet center line coincides with the neutral 
axis of the ring. Of major interest in a structure such as the one 
described herein are the forces, stresses (obtainable from the forces), 
and displacements listed in the following paragraphs, together with their 
mathematical expressions which are readily obtained in general form in 
the appendix. These expressions become considerably simplified in the 
absence of distributed loads. (For positive sign convention of the 
displacements and forces see figs. 1 and 3# respectively.) 


The eight quantities listed in this paragraph are associated with 
the boundary conditions (discussed in the next section). These 
quantities - the forces and displacements at the flexible ring, the 
shear flow at the panel edge, and the axial displacement between the 
flexible and rigid rings at the panel edge - are: 


( 1 ) Bending moment, M 

( 2 ) Shear force 


( 3 ) Axial force 


„ 1 dM 

-I® 


R d 02 d 0 


(4) Sheet shear flow 


q = 


2R^^ 


j ^ _ 

I ^ d0^ 


R' 


(5) Axial displacement 


u = 


4R3Et' 


d^ d\l „2 

-p + — r - ^ 

d 0 


df 1 d 

de R r ^2 , 


2 / 3 ^ 

dh ^ ^ ^ 1 /to ^ ^ 


de 


d 0 


R\d0 


de- 
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(6) Tangential displacement 
2R Gt 



— — 1 


dM d% „2 

, , df , 1/ d^ 


ds 3 

t ^ 

> 

L 




J d%t ^ ^ _ 


6R^t' |dj03 d05 

( 7 ) Radial displacement 


R^ 


d^h ^ d^f ^ 1 A^m dSn^ 


de-^ 


dsS 



w = - 


L , 

fd% d^ „2 

dh d^f 1 /4a d4A 

2R%t 

aea * " 

>h. 

dS ^ d0^ ^ R^dS de^ 


> + 



V. 

foi H 

d^ 


d^h j d^f. 1 l/4^m 1 

-V 

6R%:t' 

Ijie^ 

d0^ 


d0^ dS^ ^\de3 ds^y 

> 


^ - e 

EI\ R d02 dey^ 

(8) Rotational displacements 


dw ^0 
ds R 


2R3&t 


^ + 2 - R^ 

^ de^ ds5 


i£ + a ^ 


dB‘ 


d9 




— *lf + 

1_ dea 

d^M 


6R5Et' 


d£f 

ds3 


_ + 2^ + ^ 
d03 d©5 gsT 


Vw 


^ + ^ + £f + £f + + 2 


d0' 


d0 


^ de3 de5 


^(35 ‘ 



e AdM e d^M 

^\^de R dg3 


df d^m 


dB‘ 


The shear and normal stresses in the two supporting panels follow. 
The shear stress t in the sheet is given by 


T = 


1 

dM 2 

V cLf , 1 / , d^iA 


2R^t ^ 


^ d07 

' 






(6) 
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The stringer normal stress o is 
L - 


a = 


r u 

— 0 / q \” 

“V 

Jd^M ^ d M 2 



d02 

63 (jg2 qqZJ 



— 



2R3t' 


The ring tending moment, that Is, the bending moment about the 
neutral axis of the ring, is 

® - eRf - e ^ 


(T) 


Mr = M + - 

^ R d02 


d6 


( 8 ) 


For no eccentricity the moment at the sheet center line becomes the 
ring bending moment. 


BOUNDARY CONDITIONS 


Inasmuch as a general differential equation of the form of 
equation (l) applies for each of the panels comprising a given structure, 
a like number of bending -moment solutions of the form of equation (3) 
or ( 5 ) results. For the determination of the unknown coefficients 
appearing in these solutions a sufficient" niamber of boundary conditions 
must be found, one for each coefficient. Regardless of the number of 
panels, the boundary conditions involve only the eight expressions 
listed in the previous section, consisting of the forces and displacements 
which occur at the sheet center line. 

Although, for a particulaj: problem, the application of these 
expressions depends upon the structiire to be analyzed and the loading, 
all the boundary conditions may be summarized in one general statement; 

No change in displacement or forces can occur across a boundary unless 
a concentrated force is applied at the boundary, in which case the change 
in forces must equal the applied force. Any concentrated load in the 
plane of the ring can be resolved into a radial or shearing force, a 
tangential or normal force, and a moment. Then when a concentrated load 
is applied, the boundary conditions require that the difference in 
shear forces of the adjacent panels be equivalent to the applied radial 
load, the difference in axial forces be equivalent to the applied 
tangential load, and the difference in moments equal the applied moment. 

In the absence of any concentrated loads, all the forces and displacements 
must be continuous; that is, all eight expressions in one panel must 
equal the corresponding expressions in the adjacent panel. If a concen- 
trated load is applied within a panel, it is necessary in the analysis 
to consider the point of application of the load as a boxindary and, hence, 
to consider the panel as two panels, one on either side of the load. 
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Although terms associated with distributed loads appear in the 
expressions for the boundary conditions, distributed loads affect the 
boundary conditions only indirectly inasmuch as they affect the 
displacements . 

Further discussion of boundary conditions, as they apply to one of 
the n\imerical examples, is given in the following section. 

HUMERICAL EXAMPLES 

The numerical examples deal with two geometrically similar cylinders 
constructed of four panels forming a doubly symmetrical, essentially 
elliptic cross section (fig. 2), with the flexible ring subjected to 
concentrated radial, moment, and tangential loads at an intersection 
with the major axis. Cylinder 1 has a very stiff loaded ring and 
cylinder 2 a relatively flexible loaded ring. For each cylinder the 
sheet thickness is constant, there is no eccentricity of ring and sheet, 
and the radliis of the top and bottom panels is one -third the 

radius R 2 of the middle panels. The moment of inertia of each ring is 

constant, but becaxise of the change in radius, the rings change in 
relative stiffness from one section to another as Indicated by the change 
in a/b given in table II. As seen in figure it, the top panel, section 1, 
Joins the middle panel, section 2, at 0 = 60°, and the middle panel 
Joins the bottom panel, section 3> at 0 = 120°. (These dimensions were 
also used in constructing fig. 2.) 

Comparisons of the calculated distributions of bending moment, shear 
force, and axial force in the ring and of shear flow in the sheet with 
the distributions given by the elementary theory are shown in figures 5 
to 10. The necessary numerical values used in the calciilations are given" 

in tables II and IH. For ease in reading figures 5 to 10, the abscissa, 

although it actually represents distance along the perimeter, is given 
in degrees measured from the vertical axis of symmetry (as shown in 
fig. it). Thus, since the ring perimeter of section 2 (fig. it) is three 
times that of section 1 or section 3 ) whereas the angular distance for 
all three sections is the same (6o°), the distance in figures 5 to 10 

along the abscissa from 120° to 6o° (corresponding to section 2) represents 

three times the distance from l80° to 120° or from 60 ° to 0° (corre- 
sponding to sections 1 and 3). 

For the numerical examples, the labor of computation necessary to 
calculate the bending moment and other desired quantities may be shortened 
somewhat through cognizance that antisymmetrical loading produces an 
antisymmetrical moment distribution about the vertical axis and symmetrical 
loading yields a symmetrical distribution. Hence, only half the cylinder 
at the ring need be considered (fig. ^) . The procedure used in obtaining 
the numerical results is illustrated by taking as an example cylinder 2 
subjected to a concentrated radial load at 0 = 180°. The discussion 
to follow is confined to cylinder 2 so loaded. 
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Differential Equations and Solutions 
In the top panel, section 1 (fig. 4), the differential equation is 

A _ ^ A ^ ^ _ 20 ^ 900000 „ _ o 

^8 3 (je6 3 (304 "T 302 729 

The resulting auxiliary algebraic equation yields the eight complex 
roots (see table II): 

±(ci ± dj,l) = ±(2.509792 ± 0.6216651) 

+(c 2 + dgl) = ±(0.915513 ± 2.1019991) 

Since symmetrical loading is applied to the structure, the moment in 
this section is given by only the symmetric terms: 

Mgectlon 1 ~ ^2 sinh c^0 sin d^S + cosh cos d^^S + 

Cg sinh CgS sin d .26 + C 3 cosh CgS cos d 2 S (9) 

In the middle panel, section 2, the differential equation is 

^ - 58 ^ - 119 ^ - 60 — + 900000M = 0 

d0° de° d0^ d0^ 

The aiixillary algebraic equation has the eight Toots (see table II): 

±^L = ±5.534415 

±^2 = ±7.395633 

±(a ± bi) = +(2.181206 ± 4.2918821) 

The moment in section 2 is then 

^section 2 “ ^ 9 ' sl^lh + C^o' cosh + C 2 _ 3 _' sinh kg 0 + 

C 12 ' cosh k 20 + sinh a0 cos b0 + 

sinh a0 sin b0 + C 3 _^ cosh a0 sin b0 + 

C 16 cosh a 0 cos b 0 
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"but since ki 0 and , when 9 = 6o° and 9 = 120 °, are large enough 
so that their respective hyperbolic sines and cosines are of almost 
equal magnitude, it is better for computational purposes to rewrite the 
moment in section 2 as 

k-i 9 -kn 9 kpS -kp0 

^section 2 = ^96 1 + C^qS + C^^e + 0^2® + 

C-[_2 slnh a0 cos b0 + sinh a0 sin b0 + 

*^15 ^ ^16 cosh a 0 cos b 0 (lO) 

For the bottom panel, section 3 f the differential equation is the 
same as for the top panel and the axudliary algebraic equation has the 
same roots. Hence, the moment is given by 

Mgectlon 3 “ *^17' <^1® cos d^_0 + 02^0’ slnh C]_0 sin d^^0 + 

C19' cosh 02^0 sin d2^0 + C^q' cosh C2_0 cos d2^0 + 

C p-] sinh C20 cos d20 + C22 sinh 02© sin d ^9 + 

C20 cosh C20 sin dg0 + C2IJ. cosh CqB cos d^ 

However, for 0 = 120 ° and 0 = l 80 ° the hyperbolic sines and cosines, 
respectively, of C3_0 are almost identical, and it is advisable for 

computational purposes to rewrite the moment in the bottom panel in the 
form 

Cj_0 -Cj_0 

^section 3 = ^17® <^1® + C20e cos d^_0 + 

Cj_0 -C2^0 

C^^^e sin dj^0 + CgQ® sin dj^0 + 

C21 slnh C2B cos d2© + C22 sinh c ^9 sin d20 + 

C23 cosh C 20 Bin d 20 + C24 cosh C20 cos d20 ( 11 ) 


Boundary Conditions 

Equations (9) to ( 11 ) contain twen"ty unknown constants; hence, 
twenty boundary conditions are needed.. No boundary conditions are found 
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at 0=0, since continuity of all forces and displacements is already 
satisfied as a consequence of taking advantage of the symmetry of the 
structure and loading. The boundary conditions to be used in the 
calculations must be found at 0 = 6o°, 0 = 120°, and 0 = l80°. 

Eight of the boundary conditions are supplied at 0 = 60°, where, 
in the absence of any concentrated loads, all eight expressions in 
section 1 must equal the corresponding expressions in section 2. That 
is, the moment, shear force, displacements, and so forth, must be 
continuous . For example , the first condition requires continuity of the 
moments, or 


^section 1 “ ^section 2 ^ 


which can be written 


c-i n ci Jt dn It 

C 2 slnh sin -f C 2 ^ cosh cos -t- -f 


Colt doic 

Cg sinh — ^ sin + 



Cplt 

Cg cosh — ^ cos 


kj^Jt 






kg It 
1 ” 




d2* 

1 “ 

3 . 



sinh 

cosh 


art 

T 


art 

T 


14 

16 



0 


The other seven boundary conditions appearing at 0 = 6 o° are written 
in a similar manner. The same continuity of the forces and displacements 
must exist at 0 = 120°. Thus sixteen boundary conditions are provided, 
with the four remaining conditions to be found at 0 = l 80 °. 


For an applied radial load, the sum of the shearing forces on either 
side of the cut must equal the applied shearing load P. Because of 
the symmetry of the structure and loading, half of this load is carried 
by section 3 and the other half of the load is carried by the peinel 
to the left, which need not be considered. Furthermore, the shear 
flow in the sheet is 0 at 0 = l 80 °, there can be no tangent leil or 
sidewise displacement at 0 = l 8 o°, and there are no angular displace- 
ments caused by a concentrated radial load at this point. Hence, the 
four boundary conditions needed may be summed up as follows; 
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^ 180 ° - 0'5P 

‘llSoo = ® 


(12a) 

(12Td) 


^0l80° = ° 


(12c) 



= 0 


(I2d) 


Equations (l2a) to (I2d) can "be reduced "by proper sulDstltution to 


(a) 

Vd9 y n 

^ '^180° 

0.5PR1 

(^^) 


) ■ 

= -0.5PR1 


\d0 y 

' ^ 0 
180 


(c) 

/d 5 M^ 

0 

0 

= 0.5RR1 

(d) 

/'d'^M^ 
\ 7 i 

1 

= -0.5PR1 


\d0 7 

180° 



For example, equation (12a) can "be written 


^1^ y N ”'"1 ‘'"I \ 

(c^ cos d 3 _s - di Bin - C^ge (c^ cos + d^^ sin d^itj 


-CtK , 


X \ 

sin djit + cos ^1^) 

C 21 (cg cosh C2^ cos ^ sinh C2^ Bin ^2^ + 

^22(^2 ■*■ ^ ^ 2 ^ 

C23 (cg sinh CgK sin d£« + ^ cosh Cg* cos + 

C24(sinli C2« cos ~ ^ cosh Cg« sin 62^ = 0 . 5 Hl-]_ 


-CtK , 
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The solution of the twenty equations given by the boundary conditions 
yields the values of the unknown coefficients Cg, . . . C23, C2IJ. 

shown in table III. The bending moment in the loaded ring at any angle 6 
is then found by using equation (9)^ (lO), or (ll) with the appropriate 
coefficients. For example, substituting coefficients C-j^^ to C^^ 

from table III into equation (ll) for 9 = l 80 ° gives 

MiqqO = -0.11905^ X io“^pRi( 2.656MA- X lo^C -0.372983) + 

(-0.202194 X 10 ^Ri)(). 376443 X 10"5 (-0.372983) + 

^0.342801 X 10 (2. 656444 X io^}(^o. 927839)+ 

(0.216297 X lO^PR^ (0.376443 X 10"^) (0.927839) + 

(^0.497509 X 10 “^PR^^( 8 . 444906) ( 0 . 949098) + 

(■0.226239 X 10 ~%Ri)( 8 . 444906) (0.314983) + 

(0.3751^48 X l0‘^PRi)( 8 . 901257) {0.314983) + 

(0.532861 X 10 "^ERi^( 8 . 901257) (0.949098) 

= 0.161390PR1 


CONCLUDING REMARKS 

The resvilts obtained from the numerical examples agree with those 
previously obtained for circular cylinders (see, for example, references 2 
and 3) in indicating that concentrated loads applied to flexible rings 
produce stresses in the rings and shell considerably different from those 
computed from an engineering analysis (wherein the ring is treated as 
a free ring supported by the usual elementary torsion and bending shears). 

Ring flexibility is essentially Indicated by the parameter A I = — =— Ji 

V IL-^ / 

where A is everywhere less than about 200 the engineering analysis is 
adequate, but if A exceeds 200 such an analysis is inadeqtiate. 
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The main effects of flexibility are to change considerably the distri- 
bution of stress and the magnitudes of the maximum stresses. The change 
in maximim stresses is indicated in the following table, which gives the 
approximate ratios of the absolute values of the maximum stress coeffi- 
cients for cylinder 2 (A 206 near the load but 150,000 some distance 
away) to the Tnaximum stress coefficients obtained from an engineering 
analysis: 



For much larger values of A, that is, for greater flexibility, the 
ratios greater than unity would increase considerably, those less than 
unity would decrease somewhat, and the ratios of unity would remain 
unchanged. 


Langley Aeronautical Laboratory 

HatloneJ. Advisory Committee for Aeronautics 
Langley Field, Va., A\igust 1, 1951 
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APPENDIX 

DEVELOPMENT OF GENERAL DIFFERENTIAL EQUATION 


Consld.©!* "the 'two-Uay panel vith. unlfonn geome'ti'lcal properties 
(^^S* 1 )^ loaded in the plane of the middle ring. In accordance with 
the assumption that the shear stress in the sheet may vary circumferen- 
tially hut remains constant in the axial direction, an infinitesimal 
element of the sheet is subject to the stresses shown in figiure 1. 

Equilibrium in the x-dlrection requires that 


Since 


— t ds dx + — t' dx ds = 0 
os ox 


(Al) 


and 

eqxxation (Al) becomes 


Tt = q 

crt' = p 




(A2) 


Integrating equation (A2) with respect to x gives 


P = 



fl(s) 


However, since the end ring supports are free to warp out of their planes, 
p=0 at x=L, and therefore 


fl(s) 



Then 


P 




The strain relation in the x-dlrectlon gives 
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so that 


i (L - x) 

Sx Et' Ss 


Integrating equation (A3) with respect to x gives 


u = ^(Lx-^)|a + f2(8) 
Et* V 2 / Ss 


hut since, from, symmetry of the structure and loading, u = 0 at x - 0, 

f2(s) =0 


so that 




The shear strain can he expressed as 

y = ^ ^ (A5! 

Gt Ss 5x 

from which u can he eliminated and v can he found in terms of q. 
Differentiating equation (a 4) with respect to . s and suhstltuting the 
result into equation (A5) gives 

^ = -9- - -A- (lx - 
5x Gt Et* V 2 / Sg2 


Integrating with respect to x yields 

q 1 /Lx^ x^'^ S^q . ^ f s 

hut since the- end ring supports are rigid in their planes, V = 0 at 
X = L, and 


SO that 




I- (3I8C^ - x3 - 2 l3) ^ 
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This displacement in the sheet at the ring, where x = 0, hecomes 


^0 


_L_ 

Gt 


tl + 


l3 

3R^t' 


(a6) 


The relation between the strain in the sheet in the circumferential, 
direction and the deflections at the ring is 




^^0 

W~ 


w 


where by continuity w is the radial deflection of both ring and. sheet. 
Since the strain Cg in the sheet is the strain at the outermost fibers 

of the ring, 

6 = 

® E ' El 


Thus the radial deflection can be expressed as 


w 




(AT) 


Differentiating equation (a 6) and substituting the result into 
equation (AT) yields 

- ir- - Gt d0 ^3 


(A8) 


The relation between the ring bending moment and the radial 
deflection is given by (see, for example, reference 11) 


Mr = - 



(A9) 


Performing the indicated differentiation of equation (A8) and substituting 
the result into equation (A9) eliminates w and gives 


Mr = - 


El 

(R - e)2 




+ 



L d^ ^ d£A 

Gt^g3 flRy 3R^t'\^^ d©y 


(AlO) 
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It is now convenient to find the relation between the shear flow 
in the sheet and the moment at the sheet center line, as well as to 
express the moment at the ring center line in terms of the moment at 
the sheet center line. Figure 3 shows an infinitesimal segment of the 
loaded ring and the positive directions of the forces. For the first 
of the three equations of equilibrium, taking moments about the origin 
gives 

-dM + R dE + dS + Rm cl0 + 2 r2<i dS = 0 


or 


q = 


1 fm 

2R^^ 




(All) 


Summation of the tangential forces (neglecting terms of higher order) 
gives 


V <ae - Rh dR - 2Rq dfi - dH » 0 


or 


V = ^ + Rh + 2Rq (A12) 

Summation of the forces in the radial direction gives 

-Hdfi+Rfde-cW=0 


or 


H = - ^ + Rf (AI 3 ) 

Substituting the value of from equation (A12) into equation (All) 

CL“ 

yields the shear force at 0; 

(AlW 

Differentiating equation (Al4) and substituting the result into 
equation (AI 3 ) gives the value of the normal force at the sheet center 
line at 0 : 
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Differentiating equation (A15) and subtituting the result into 

equation (All) gives the expression for the shear flow in the sheet at 6: 


q = 




h + 


as 



(Al6) 


The bending moment about the ring' center line can be seen to be 
composed of the moment M and the product of the axial force H and 
the eccentricity e, or 


Mr = M - eH 
= M 

R de2 


eRf - 



(A17) 


Differentiation of equations (A16) and (Al?) and substitution of 
the result into equation (AlO) leads to the desired differential equation: 


A ^ (2 . 3B) A ^ 6e5A^ A . 


d0 


8 


d0 


B - 2 - ^ ^ 

R V R r£ 




dSl 


d©' 


+ 6a(1 - - + — )M = F(0) 


r2> 


(A18) 


where 


F(0) = R'^ 


^ + (1 . 3B)^ . 3B 

d© 


d©5 


.d©3 


Id^f 


>d^f 


^4 + (1 . 38)^4 - 3(b - Ssfy 5% + M(i . £ + £!)fi 

^ aa** V rV ae® ® \ e7 



, 4. (2 - 3b) 3^ + (i - 6 b * - 3 

' d©^ V R‘=^/d©3 


B - 


2eA| 


1 - i + SL 


I 


dm I 

d©; 


Equation (Al8) is the general form of the differential equation for the 
bending moment at the sheet center line for_any of the panels comprising 
the structure. The bending moment in the ring for any such panel is 
given by equation (AIT). _ _ 
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TABLE I.- BOOTS OF AUULIAHI ATOBHRAIC EqQfiTIQV WHEN e - 0 FOR 


SEVERAL CCMBINAnONS OF A ABB B 

t Il3 ' 



eo 

60 

150 

350 

Root 

STisbols 


±1.250205 

+0.956180 

±0.743842 

±0.571042 

±tl 

10^ 

±7.7J^5793 

±13.4 i 64 o 4 

±21,213203 

±32.403439 

+kg 


+(0.612234 + 1.4679881) 

±(0.47555^^ ± 1.2973971) 

+(0.371310 ± 1.1893351) 

+(0.290865 + 1.1131441) 

±(a + 1)1) 


±2.027^36 

±1.613689 

±1.326187 

±1.092601 

I’^l 

10^ 

±7.744229 

±13.416370 

±21.213202 

±32.403439 

±kg 


±(0.972538 f 1.9969011) 

+(0.798821 ± l.71‘^525l) 

±(0.661073 + 1.5219621) 

±(0.547534 + 1.3740901) 

+{a ± tl) 


±3.157270 

±2.523294 

±2.117014 

±1.795392 

+4± 

10^ 

+7.728363 

±13.416028 

+21.213188 

±32.403439 

+^2 


±(1.447531 ± 2.8183901) 

±(1.23639^ ± 2.3889511) 

±(1.052130 + 2.O85072I) 

+(0.896469 ± 1.8457881) 

+(a ± M) 


±5.016203 

±3.835266 

±3.230542 

±2.769566 

% 

10^ 

±7.539604 

±13.^12609 

±21.213049 

±32.403431 

+^2 


±(2.058771 ± 4,0301991) 

+(1.837706 ± 3.4177221) 

+(1-595017 + 2.9595621) 

+(1.379274 + 2.5939991) 

+(a ± bl) 
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TABLE II.- NUMERICAL VAHJES USED IN COMPUTATIONS 


A - • B =1^ 

Il 3 ' ® GtL^ 
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TABLE III.- VALUES OF COEFFIdEHTB^ FOR NUMERICAL EXAMPIES 


MM 

Cylinder 1 

Cylinder 2 


Load 

Load 

■M 

Radial 

Moment 

Tangential ^ 

Radial 

1 

Moment 

Tangential 



0.147507 

0.249767 


-0.860733-03 

0.119301-03 

^1 




0.679031-01 



0.949137-04 



Lq 







.275041-04 

-.157406-03 


.292731-03 

-.371518-04 


.125895-05 

-.284700-03 











-.136156 

.754809 1 


-.148018-01 

.221238-02 j 


-.525229 

-.197121-02 











-.165618 

.236338-03 


-.521539-02 

.484959-03 





.282018-01 



-.354900-02 








°9 

.171565-01 

.188166 

.20246? 

-.657213-08 

.14 o 829 -o 6 

-.161423-07 

=10 

.303490-01 

-.121785 

-.111398 

-.534558-01 

.270768 

-.697692-01 

°11 

-.264559-10 

.116591-10 

-.116591-10 

.397640-09 

-.161540-08 

.710774-09 

°12 

.889303 

-.508263-01 

.851830 

-.479259 

-.380612+01 

.638771 

°13 

.329176-01 

-.296012 

.472934 

-.122989-01 

-.361474-01 

.485865-02 

Ci4 

-.490646 

.373871-01 

.290987 

.829422-03 

-.491296-01 

.702908-02 

°15 i 

.192013 

-.184225 

-.160403 

-.117551-02 

. 471184-01 

-.668234-02 

Ci6 

.720834-01 

.274929 

-.400188-01 

.131739-01 

.392728-01 

-.527300-02 

°17 

-.119460 

-.164389-01 

.174960 

-.119054-05 

.778086-04 

-.192396-04 

C18 

.315057 

.363960 

.224540 

-.202194+01 

-.711069+01 

.872951 

0 

H 

-.898224-08. 

.174922-09 

.301271-07 

.342801-04 

.8o4877-o4 

.168139-04 

°20 ■ 

-.711275-02 

-.259599 

-.364881 

.216297+01 

.882602+01 

-.112279+01 

Cgi 

-.640348-01 

.207820-02 

.128917+01 

-.497509-01 

-.142194 

.194136-01 

°22 

-.531867 

.131177-01 

.506058-01 

-.226239-01 

-.709217-01 

.370731-02 

C23 

-.254629 

-.304436 

-.328941 

.375148-01 

.124468 

-.104012-01 

C24 

.194401-01 
i 

-.2lte335“01 

1 

-.613555-03 

.532861-01 

.167720 

~. 2 kl 62 hJ^l 

\ 


^he following conyention is used to indicate a\iltiplication factors: 



tOl * X lOj +02 » X 10 f -01 “ X 10 -02 » X 10 } nTvi BQ forth. 
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Figure I.- Two-bay panel showing positive directions of displacements 
and sheet stresses. 



Figure 2.- Two- bay cylinder , composed of four clrculor panels, with doubly 
symmetrical noncircular rings, used for numerical examples. 
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Figure 3.- Free-body diagram of ring section showing positive 
directions of forces. 





Figure 4.- Cross- sectionai view of half of loaded ring used for 
numerical examples. 
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— Elementary theory 




Figure 3.- Bending -moment and shear-flow distributions produced by 
concentrated radial load P at 6 = 180". 
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— Elementary theory 
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(b) Sheet shear flow. 


Figure 7.- Bending-moment and shear- flow distributions produced 
by concentrated moment load Mg at » I8O”. 



(b) Ring axial force. 


Figure 8.- Ring shear- and axial- force distributions produced by 
concentrated moment load Mq at ■ 180“ 
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— Elementary theory 

Cylinder I I _ 

Cylinder 2 J 
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(a) Ring bending moment. 



I 
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I 


100 80 
S, deg 

Cb) Sheet shear flow. 


Figure 9.- Bending -moment and shear-flow distributions produced 

by concentrated tangential load T at 5 ■ 180® 
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— Elementary theory 

Cylinder I 1 

— Cylinder 2 J 




Figure 10.- Ring shear- and axial-force distribLitions produced by 
concentrated tangential load T at fl»l80®. 
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